A Conserved variable in the perturbed hydrodynamic world model 



Jai-chan Hwang 

Department of Astronomy and Atmospheric Sciences, Kyungpook National University, Taegu, Korea 

(February 1, 2008) 

We introduce a scalar-type perturbation variable $ which is conserved in the large-scale limit con- 
sidering general sign of three-space curvature (K), the cosmological constant (A), and time varying 
equation of state. In a pressureless medium <!> is exactly conserved in all scales. 

PACS numbers: 98.80.Hw, 98.80.-k 
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1. Introduction: Relativistic cosmological perturba- 
tions with hydrodynamic energy-momentum tensor was 
originally studied by Lifshitz in 1946 Q based on the 
synchronous gauge. More convenient analyses based on 
better suited gauge conditions were made by Harrison 
in 1967 U and Nariai in 1969 j|. Although the vari- 
ables used by Harrison and Nariai are free of gauge 
degree of freedom (thus, equivalent to corresponding 
gauge- invariant variables), these gauge conditions be- 
came widely known by a seminal paper by Bardeen in 
1980 0]. We believe these are brighter side of the history 
concerning gauge condition in the cosmological pertur- 
bation: Lifshitz carefully traced the remaining gauge so- 
lutions, Harrison and Nariai in fact hit the correct gauge 
conditions for handling aspects of hydrodynamic pertur- 
bations [see below Eq. (|)], and Bardeen showed the 
gauge invariance of the variables used by Harrison and 
Nariai and demonstrated the diversity of gauge condi- 
tions and gauge-invariant way of handling them. On the 
other side, there exist persistent algebraic errors in the 
literature which are often claimed to be due to wrong 
gauge conditions These errors probably gave some 
researchers an (inappropriate) impression that the field 
is 'plagued with gauge problems'. In any case, there ex- 
ist many gauge conditions waiting to be employed with 
possible advantages in exploring certain aspects of prob- 
lems. In 1988 Bardeen || made a practical suggestion 
concerning the gauge condition which allows maximal 
use of the various different gauge conditions depending 
on problems. In || we elaborated the suggestion and 
recently termed our approach a gauge-ready mothod; see 
Sec. 2. 

We may not need to emphasize the importance of con- 
served quantities in physical processes. In cosmological 
perturbations the conserved variable can provide easy 
connection between the final results and the initial condi- 
tions. Aspects of conserved perturbation variables were 
discussed in pQL|l4]|. In this paper, based on the gauge- 
ready method, we will derive a scalar-type perturbation 
variable which is conserved in the large-scale limit in- 
dependently of changing background world model with 
general K, A, and the perfect fluid equation of state. 

Besides the scalar-type perturbation we also have 
the vector-type (rotation) and the tensor-type (gravita- 



tional wave) perturbations which evolve independently 
to the linear order in our simple background world 
model. We also have conserved quantities for these 
additional perturbations: ignoring anisotropic stresses, 
the angular-momentum of the rotation is generally con- 
served, whereas the non-transient solution of the gravi- 
ational wave is conserved in the super-horizon scale in 
near flat case. These two conservation properties were 
already noticed in pj, and recently, we have elaborated 
these conservation properties in some general situations 
p"5| . In the following we concentrate on the scalar- type 
perturbation with the hydrodynamic energy-momentum 
tensor. 

Although redundant, in order to make this paper self- 
contained, in the Appendix we present the complete set 
of perturbed equations based on Einstein equations. Our 
new result is the conserved variable in Eq. (|9|) with the 



equation and the large-scale solution in Eqs. ( |13| , |ll[ ). 

2. Notations and strategy: We consider the most gen- 
eral scalar-type perturbations in the spatially homoge- 
neous and isotropic world model. Our notations for the 
metric and the energy-momentum tensor are: 
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where = rj, and an overbar indicates a background 
order quantity and will be ignored unless necessary. Spa- 
tial indices (a, /?,...) and V Q are based on g^ which is 
the three-space metric of the homogeneous and isotropic 
space. (3 and 7 always appear in a spatially gauge- 
invariant combination \ = a(f3 + a'y)] an overdot denotes 
the time derivative based on t with dt = adrj. Using Xi 
all the perturbed metric and energy-momentum tensor 
variables in Eqs. (]lj,@) are spatially gauge- invariant. 

Perturbed order variables can be expanded in eigen- 
functions of the Laplace-Beltrami operator (A based on 

dad) w ith eigenvalues —k 2 where k is a comoving wave 
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number 0; to the linear order each eigenmode decouples, 
and without any confusion we can assume variable in ei- 
ther the configuration space or the phase space. For the 
flat (K = 0) and the hyperbolic (K = — 1) backgrounds 
k 2 takes continuous value with k 2 > 0, whereas in the 
spherical (K = +1) background we have k 2 = n 2 — K 
(n = 1,2,3,.. .) 0f|m|/§. Situation in the hyperbolic 
background may deserve special attention. Probably be- 
cause any square integrable function can be expanded 
using harmonics with k 2 > 1 (subcurvature) modes only, 
most of the cosmology literature ignored < k 2 < 1 (su- 
percurvature) modes, and gave a wrong impression that 
the supercurvature modes do not exist [Qj^]: a state of 
affairs was well summarized in JlTj . In the spherical back- 
ground, the lack of physically relevant perturbations for 
the lowest two harmonics n = 1 (k 2 = 0) and n = 2 
(k 2 — 3K — 0) was pointed out in jlj. Later, we will 
use vanishing c 2 k 2 term as the large-scale limit. Thus 
the results in such a limit will be relevant in all scales for 
general K in the pressureless medium (c 2 — 0), and in the 
large-scale limit (k 2 — > 0) for the flat and the hyperbolic 
situations in a medium with c 2 ~ 1(= c 2 ). 

Complete sets of equations for the background and the 
perturbed orders are presented in the Appendix. Equa- 
tions ( ^9[p5| ) are written in a gauge-ready form. In this 
form we still have a right to choose a temporal gauge 
condition, fl. Equations are designed so that impos- 
ing the gauge condition is as simple as the following: 
the synchronous gauge chooses a = 0, the comoving 
gauge chooses v/k = 0, the zero-shear gauge chooses 
X = 0, the uniform-curvature gauge chooses ip = 0, the 
uniform-expansion gauge chooses k = 0, the uniform- 
density gauge chooses 6/j, = 0, etc. Except for the syn- 
chronous gauge, since any of the other gauge conditions 
completely fixes the temporal gauge degree of freedom 
[see Eq. (g7|)], any variable in such a gauge condition 
is equivalent to a unique gauge-invariant combination of 
the variable concerned and the variable used in the gauge 
condition. Examples of some useful gauge-invariant com- 
binations can be constructed using Eq. (p7|): 
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In this way, we can flexibly use the gauge degree of free- 
dom as an advantage in handling problems PJy]. 

3. A conserved variable: From Eqs. (|lf 
can derive JI§| ] 
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In the large-scale limit, super-sound- horizon scale where 
we ignore c 2 k 2 /a 2 term ]19| , and ignoring stresses (e and 
a), Eq. (^) has an exact solution valid for general K, A 
and time- varying equation of state p(fi) , pOp] 



^ x (x,t)=C(x)4^G 
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dt + d(x 
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where C(x) and d(x) are coefficients of the growing and 
decaying solutions, respectively. ip x most closely resem- 
bles the behavior of perturbed Newtonian potential 0. 
The variables most closely resembling the Newtonian be- 
haviors of the density and velocity perturbations are 5 V 
and Uy, respectively From Eqs. (^0 21), and Eqs. 

( |T9| , pl| , p2| ) , we can derive, respectively: 
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Solutions for 5 V and v x follow from these equations. We 
can similarly derive a solution for tp v using 
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where we used Eq 
er, we can show 

$ = 
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Thus, using the large-scale solution in Eq. (^|) we have 

d>(x,t)=C(x), (11) 

where the decaying solution has vanished. Therefore, $ is 
generaly conserved in the super- sound-horizon scale con- 
sidering general K, A and time-varying p(fi). On the 
other hand, from Eqs. (|],[l0|), thus ignoring stresses, we 
have 
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Thus, in a pressureless case <&(x, t) = G(x) is an exact 
solution valid in general scale; this was noticed in p3| . 
Combining Eqs. ( |i7]| , |l2| ) we have a closed form equation 
for $ 
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which is valid for c 2 ^ 0. Thus, for vanishing c 2 k 2 /c 
term [[l9| we have a general solution 



$(x,*) = C(x)+d(x) 
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which includes c 2 = limit. We can show easily that 
d term in Eq. (|l4|) is higher-order in the large-scale ex- 
pansion compared with d term in Eq. (^): by comparing 
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the two solutions of Eqs. (|,|T|) in Eq. (|T|) we can show 
d = -k 2 d — Ad. Therefore, the solution in Eq. ( pi] ) is 
valid in the large-scale (super-sound-horizon scale) with 
vanishing dominating decaying solution. 

4- Other conservation variables: $ differs from the 
well known conserved variable C in [fiOp|]. In £ was 
introduced in the flat background and in that background 
it is the same as 
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whereas, in the flat background we have <E> = tp v Q . It is 
interesting to note that (ps is also conserved in the large- 
scale limit considering general K jjj: from Eqs. flPlllHl) , 
we can derive 



<A5 = $ 



1 



12irG((i + p) a 



k 2 



(16) 



According to the solutions in Eqs. (ml l[) , or Eq. (|bj), we 
can see that for vanishing k 2 order term (ps is conserved 
considering general K. However, we also see that due to 
the second term the conservation property of ipg breaks 
down near and inside horizon (k/[aH\ > 1) even for K — 
0; whereas, $ is conserved independently of the horizon 
crossing in the matter dominated era. We can also show 
the conservation property of tp K : from Eqs. ( |27| , |20| ) we 
can show 
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Thus, the conservation property of (p K breaks down for 
K ^ or near and inside horizon. In the K = sit- 
uation, (p in many different gauge conditions shows the 
conserved behavior in the large-scale limit: for the ideal 
fluid see Eqs. (41,73) in §§, and E Q- (34,35) in §4|; for 
the scalar field see Eqs. (92) in p5| ; and for the general- 
ized gravity see Sec. VI in EG]. Conservation properties 
of (p in various gauge conditions were also discussed in 
p3| , PH where the arguments were based on the first or- 
der equations of the type in Eq. (|l2]). 

5. Discussions: We would like to emphasize again that 
the conservation property of <E> is valid in the limit of van- 
ishing c 2 k 2 /a 2 term. Thus, in the pressureless medium 
(c 2 = 0) it applies in all scales for general K, and in 
the medium with dominant pressure (c 2 ~ 1) it applies 
in the large-scale limit (k 2 — > 0) for the flat and the hy- 
perbolic situation. As long as these conditions are met, 
$ is conserved independently of time varying equation 
of state p(n). Since we anticipate time varying equa- 
tion of states during equal-time from the radiation domi- 
nated era (p = ^/i) to the matter dominated era (p — 0), 
and during (p)reheating period from the acceleration era 
(p < — to the radiation era, and since the obser- 
vationally relevant scales stayed in the large-scale limit 



during the transitions, $ is a practically important quan- 
tity in tracing the evolution of scalar-type perturbation 
from the early universe till recent era before the nonlinear 
evolution takes over. 

We would like to conclude by remarks on the history 
about the variables in Eq. (^): ip x and S v in their gauge- 
invariant forms became widely known by Bardeen's work 
in 1980 (these are &h and e m in 0). However, ip x is 
the same as <p in the zero-shear gauge [which fixes \ = 
as the temporal gauge condition, see Eq. (||)] which was 
first used by Harrison in 1967, and S v is the same as S 
in the comoving gauge (which fixes v/k = 0) which was 
first used by Nariai in 1969 ||. For K = 0, the variable 
ip v is widely recognized as a conserved variable in the 
literature. Up to our knowledge it was first introduced 
as ip in the comoving gauge by Lyth in 1985 jnj, and later 
was used in the context of the scalar field and generalized 
gravity as the large-scale conserved variable pTJ]. 

We thank Peter Dunsby, Hyerim Noh and Winfricd 
Zimdahl for many useful discussions. We wish to thank 
P. Dunsby for encouraging us to find the conservation 
variable for general K and for hospitality during our visit 
University of Cape Town where the result was found. 

Appendix 

In the following we derive Einstein equations based 
on Eqs. (0^). To the background order, from Tq and 
T£ — 3Tq components of Einstein equations and T^. b = 0, 
respectively, we can derive: 
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H z = — n- — + -, H = -4ttG((i+ P ) + —, 
(i=-3H(fi + p), (18) 

where H = —. To the perturbed order we can derive: 
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k 2 -3K TT A „ . 
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K x = 12ttG (h + p) -v, 
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k + 2Hk =(K~ 3H^j a + 4irG (l + 3c 2 s ) fj,S 

+12nGe, (23) 
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where we have introduced 

6 P (k,t)=c 2 s (t)S^(k,t) 
6= — , w(t) = -, 



+ e(k,t), 

cm = i. 

A* 



(26) 



In Eq. (19) we introduced a variable k which is the 
perturbed part of the trace of extrinsic curvature; sim- 
ilarly, for meanings of the other perturbation variables, 
see Sec. 2.1 of @. Equations (|0f§3|) follow from T °, T°, 
T B ~ \ 5 a T ^ and T S ~ t q components of Einstein equa- 





and T? , — 0, respectively. This set of equations were 



tions, respectively; and Eqs. (|ij25) follow from Tq 

This set of equations 
(41-47) of @, see also Eqs. 

£ a , we have 



originally derived in Eqs 
(22-28) in §. 

Under the gauge transformation x a 
(see Sec. 2.2 in M): 
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[1] E. M. Lifshitz, J. Phys. (USSR) 10, 116 (1946); E. M. Lif- 
shitz and I. M. Khalatnikov, Adv. Phys. 12, 185 (1963). 

[2] E. R. Harrison, Rev. Mod. Phys. 39, 862 (1967). 

[3] H. Nariai, Prog. Theor. Phys. 41, 686 (1969). 

[4] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[5] We found a common error in several textbooks Q based 
on the synchronous gauge: Eq. (15.10.57) in Wein- 
berg Eq. (10.118) in Peebles, Eq. (11.5.2) in Coles and 
Lucchin, Eq. (8.52) in Moss, and problem 6.10 in Pad- 
manabhan are in errors. Notice that these are basically 
algebraic errors partly due to unnecessarily complicated 
algebra in the synchronous gauge, see 

[6] S. Weinberg, Gravitation and Cosmology (Wiley, New 
York, 1972); P. J. E. Peebles, Principles of Physical 
Cosmology (Princeton Univ. Press, Princeton, 1993); P. 
Coles and F. Lucchin, Cosmology: The Origin and Evo- 
lution of Cosmic Structure (John Wiley and Sons, New 
York, 1995); I. G. Moss, Quantum Theory, Black Holes 
and Inflation (John Wiley & Sons, New York, 1996); 
T. Padmanabhan, Astrophysics and Cosmology through 
Problems (Cambridge Univ. Press, Cam bridge, 1996). 

[7] J. Hwang and H. Noh, |astro-ph/9907063i Gen. Rel. Grav. 
in press (1999). 

[8] J. M. Bardeen in Particle Physics and Cosmology, edited 
by L. Fang and A. Zee (London, Gordon and Breach, 
1988), 1. 

[9] J. Hwang, Astrophys. J., 375 (1991) 443. 
[10] J. M. Bardeen, P. J. Steinhardt and M. S. Turner, Phys. 

Rev. D 28, 679 (1983). 
[11] D. H. Lyth, Phys. Rev. D 31, 1792 (1985). D. H. Lyth 

and M. Mukherjee, Phys. Rev. D 38, 485 (1988). 



[12] 

[13] 
[14] 



[15] 
[16] 



[17] 
[18] 



[19] 



[20] 
[21] 



[22] 

[23] 
[24] 
[25] 
[26] 
[27] 



J. Hwang, Astrophys. J. 380, 307 (1991); Gen. Rel. Grav. 
26, 299 (1994); J. Hwang and J. J. Hyun, ibid. 420, 512 

(1994) . 

P. K. S. Dunsby and M. Bruni, Int. J. Mod. Phys D 3, 
443 (1993). 

W. Zimdahl, Class. Quant. Grav. 14, 2563 (1997); Phys. 
Rev. D 56, 3357 (1997); W. Zimdahl and D. Pavon, ibid. 
57, 7607 (1998). 

J. Hwang, Gen. Rel. Grav. 30, 545 (1998). 
U. H. Gerlarch and U. K. Sengupta, Phys. Rev. D 18, 
1773 (1978); J. J. Halliwell and S. W. Hawking, ibid. 31, 
1777 (1985). 

D. H. Lyth and A. Woszczyna, Phys. Rev. D 52, 3338 

(1995) . 

In order to derive the equation for <p x using Eqs. (jH^- 
|2e| ) we simply set y — in the equations, and then each 
variable is equivalent to a gauge invariant combination 
with subindex \. Remember that ip x is the same as if 
in the zero-shear gauge which sets y = as the gauge 
condition. 

We would like to comment on our large-scale limit used 
in Eqs. (|4|l3j). Consider an equation of a form 



= 0. 



In terms of ip = \J Z/a<j> we have (/ = 1 



By ignoring c a k term in either equation we have 



i oc constant, / dt/Z. 



(28) 



(29) 



(30) 



Thus, when we ignore c 2 s k 2 /a 2 term in Eq. ( |28[ ) as the 

large-scale limit, we assume the second term in Eq. ( psj ) 

is negligible compared with the third term. 

J. Hwang and E. T. Vishniac, Astrophys. J. 353, 1 

(1990). 

After completion of a draft P. Dunsby informed us that $ 
was also introduced in |l| in the covariant formulation, 
see Eq. (16) in [ff3. We notice that authors of |l^] took 
k 2 — > 1 as the large-scale limit in the hyperbolic case, 
thus leading to different conclusions from ours: their con- 
clusion was that, except for the pressureless case, the 
conservation property of 5> is valid only in near flat situ- 
ation. 

Since £ is gauge invariant we can evaluate Eq. ([l5|) in the 

comoving gauge, thus show £ = ip v + 8 v /(3 + 3w) tp v . 

J. Hwang, Astrophys. J. 415, 486 (1993). 

J. Hwang, Astrophys. J. 427, 533 (1994). 

J. Hwang, Astrophys. J. 427, 542 (1994). 

J. Hwang and H. Noh, Phys. Rev. D 54, 1460 (1996). 

M. Sasaki, Prog. Theor. Phys. 76, 1036 (1986); V. F. 

Mukhanov, Soviet Phys. JETP 68, 1297 (1988); D. H. 

Lyth and E. D. Stewart, Phys. Lett. B 274, 168 (1992); 

J. Hwang, Phys. Rev. D 48, 3544 (1993); ibid. 53, 762 

(1996). 



4 



